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Abstract 

As direct and indirect dark matter detection experiments continue to place stringent con- 
straints on WIMP masses and couplings, it becomes imperative to expand the scope of the 
I I search for particle dark matter by looking in new and exotic places. One such place may be 
,S^ the core of active galactic nuclei where the density of dark matter is expected to be extremely 
C^ high. Recently, several groups have explored the possibility of observing signals of dark matter 
^ from its interactions with the high-energy jets emanating from these galaxies. In this work, we 
flj build upon these analyses by including the other components of the WIMP-induced gamma 
ri ray spectrum of active galactic nuclei; namely, («) the continuum from WIMP annihilation into 
' ' light standard model states which subsequently radiate and/or decay into photons and {it) the 
direct (loop-induced) decay into photons. We work in the context of models of universal extra 
dimensions (in particular, a model with two extra dimensions) and compute all three compo- 
nents of the gamma ray spectrum and compare with current data. We find that the model 
ivj with two extra dimensions exhibits several interesting features which may be observable with 
^^ the Fermi gamma ray telescope. We also show that, in conjunction with other measurements, 
QQ the gamma ray spectrum from AGN can be an invaluable tool for restricting WIMP parameter 
1* space. 
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1 Introduction 

The search for dark matter has never experienced such exciting times. Both direct and indirect 
detection experiments are starting to probe (and constrain) what many beheve to be the preferred 
region of parameter space for particle dark matter (namely, electroweak-size couplings and masses). 
There have also been several anomalies in the data from experiments on both sides that are con- 
sistent with signatures expected from particle dark matter, but may ultimately be explained away 
as either astrophysical in origin or as experimental errors. 

While these searches continue (and the debates over possible "signals" rage), it is important to make 
sure we are uncovering every stone in the search for dark matter by considering more exotic scenarios 
and/or looking in different locations in the universe. Recently, several groups have investigated the 
possibility of detecting dark matter around active galactic nuclei (AGN) using the high-energy 
jets which emanate from these objects [H [3 HI [S] . AGN are thought to possess the densest 
distributions of dark matter in the universe which makes the probability of interactions between 
the jet particles (thought to be either electrons or protons) and dark matter particles in the halo 
surrounding the AGN core non-negligible. The basic idea is that collisions of jet particles with 
WIMPs causes an "up-scattering" into a heavy, charged non-standard model particle. This particle 
quickly decays back to the original jet and WIMP particles, but, in the process, a photon can be 
emitted from the initial, final or intermediate state charged particles. The flux of gamma rays 
produced in this way has been shown to possess several interesting spectral features (including a 
sharp cutoff which depends on the mass splitting between the WIMP and the up-scattered particle) 
that would discriminate it from any astrophysical backgrounds. 

With the expected high densities of WIMPs around the cores of AGN, the probabilities for annihi- 
lations of WIMPs could also be enhanced. These annihilations can give rise to gamma rays in two 
ways. The first is annihilation into lighter standard model (SM) states which then radiate photons 
or hadronize/decay into them (e.g., through neutral pion decay). The result is a continuum which, 
generally, has a kinematic cutoff at the WIMP mass, but is otherwise featureless. The other possible 
contribution to the gamma ray spectrum from WIMP annihilations comes from the loop-level pro- 
cesses which directly produce 'y + X final states (where X can be another photon, a massive vector 
gauge boson or a scalar depending on the spin of the WIMP [6j [7l [3 El HO]). Since the WIMPs 
are essentially at rest, the result of these annihilations arc mono-energetic photons produced with 
energy: 



where Mom is the mass of the WIMP and Mx is the mass of the particle produced in association 
with the photon. Naively, since these processes occur at loop-level, one would assume that they are 
suppressed compared to the tree-level continuum annihilations. However, it has been demonstrated 
for several models that the line-emission can compete with (and sometimes even dominate over) 
the continuum. The focus of most of these types of studies has been either the center of the Milky 
Way galaxy or its orbitting satellite (dwarf) galaxies. The characteristics of continuum and line 
emission (and their relative sizes), however, is strictly determined by the particle physics involved 
(and not the astrophysics). Thus, the gamma ray spectrum of other astrophysical objects should 
exhibit similar features. In any case, observation of line emissions would provide a smoking gun 
for indirect dark matter detection since astrophysical processes are incapable of producing such 



features. 

In this paper, we compute the complete WIMP-induced gamma ray flux for an AGN (including 
the AGN jet- WIMP interaction component along with the contributions from annihilations). Since 
previous studies focussed solely on the gamma ray spectrum coming from the interaction between 
AGN jets and WIMPs, this marks the first time that the complete spectrum has been constructed. 
We choose, as a particular example, to work in the context of the six-dimensional "chiral square" 
model [TT]. This choice is only for illustrative purposes and is due, in part, because the chiral 
square model has been shown to produce several distinguishable line emisssions from annihilations 
(a generic situation dubbed "the WIMP forest" ) |5] . 

The rest of the paper is organized in the following manner. First, in Section [2] we give a brief 
overview of the 6-d chiral square model and discuss constraints on the mass of the WIMP candidate 
in this model from relic density considerations. In Section |3J we compute the gamma ray flux 
from WIMP interactions with AGN jets. This calculation involves several pieces (which can be 
categorized as coming from astrophysics or from particle physics) and we will discuss each piece in 
some detail. Sectionfflcontains a discussion on the gamma ray flux coming from WIMP annihilations 
around AGN. In Section [5J we combine the various contributions to the AGN gamma ray flux and 
compare with current observations fro the Fermi gamma ray telescope. In this section, we also 
perform a numerical fit of the model parameters using the gamma ray flux from an AGN (as 
observed by the Fermi Gamma Ray telescope) along with the observed relic density and the most 
recent constraints from direct detection experiments. Section [6] discusses the implications of the fits 
and illustrates the usefuUness of AGN gamma ray spectra not only as a way to look for signals of 
dark matter, but as a tool to constrain its parameter space. Finally, in Section!?! we conclude. 



2 The Model 

The "chiral square" model is a model of universal extra dimensions (UEDs) in which all SM fields 
are allowed to propagate in two extra dimensions. The exact details of the model can be found in 
Refs. [m 1121 Uni IHj . Here, we only point out features that are of interest to our study. 

The two extra-dimensional coordinates can be represented by a pair of points {x^,x^) living in a 
square region with sides L. The extra dimensions are orbifolded in such a way that adjacent sides 
of the square are identified with each other, 

(y,0) = (0,y) , {y,L) = {L,y). (2) 

The folding leaves the two corners of the square which lie along the fold (at (0,0) and (L, L)) 
invariant, and identifies the remaining two corners (at (0,L) and {L,0)) as the same point. 

The Kaluza-Klein (KK) modes of the SM fields are lablelled by a pair of integers {j, k) which 
satisfy: 

fc > , J > 1 - 4,0 • (3) 



The masses of the KK modes are generally given by 



M-=Ml + ^^'^^ (4) 



where Mq is the mass of the "zero-mode" field which we identify with the SM field. The KK modes 
of fermions are Dirac particles, while the gauge fields decompose into 4-d vectors V^ and two 4-d 
scalars, V^ and V^. One linear combination of these scalars is eaten, level-by- level, by the vector 
KK modes to provide their longitudinal degrees of freedom. The other linear combination of V^ 

fi (i k) 

and V are physical gauge adjoint scalars which we denote as Vj^' . 

As is the case in the much more studied 5-d UED model, the chiral square model contains a "lightest 
Kaluza-Klein particle" (or LKP) which is forbidden (at tree-level) from decaying due to remnants of 
extra-dimensional spacetime symmetries. Thus, the LKP (which is a (1,0) excitation in the notation 
introduced above) provides an excellent candidate for cold dark matter. The exact identity of the 
LKP depends on boundary terms, but we will follow the usual reasoning and assume that these 
terms are such that colored and charged KK modes are heavier than neutral modes. Following 
these (and other) assumptions, the lightest (1,0) mode is expected to be the scalar partner of 
the hypercharge gauge boson, B\j' (which we will abbreviate as Bh)- In general, electroweak 
symmetry breaking mixes Bh with its SU{2) counterpart, Wj^' . However, the mixing angle is 
typically small and the LKP (to a very good approximation) is pure Bh- This means that the 
LKP coupling to fermions and the SM Higgs is controlled exclusively by the C/(l) gauge coupling 
gi and the hypercharge of the matter field. Of particular interest to the work performed here is 
the coupling between the WIMP, an electrically-charged SM fermion and its first KK counterpart. 
The effective Lagrangian for this interaction can be written as: 

AC = 51 [^e{YlPl + YrPr),P, + MYlPr + YRPL)ibE] Bh , (5) 

where Yl{Yr) is the hypercharge quantum number for a left-handed (right-handed) SM particle, 
^e is the SM fermion field and ipE is its KK counterpart. 

Since the LKP is a real scalar, it has mass-supppressed annihilations into light fermions and, thus, 
predominantly annihilates into pairs of heavy particles such as massive SM gauge bosons and, if 
heavy enough, top quarks. These annihilations are dominated by s-channel exchange of the SM 
Higgs boson and, as a consequence, its thermal relic density is very sensitive to the mass of the 
Higgs Mh (as well as its own mass Mb) IE]. In particular, assuming the recent discovery at the 
LHC of a spin-0 boson with mass ~ 125 GeV is indeed the SM Higgs boson, the mass of Bh is 
tightly constrained to lie in the range 190 GeV < Mb < 215 GeV. Note that this is assuming that 
Bh is solely responsible for the current density of DM. It could be that the LKP is not a thermal 
relic or that there are other quasi-stable particles which contribute to the relic density. In these 
cases, the mass of Bh may not be so tightly constrained. We will investigate the constraints on 
the WIMP mass (and other parameters) further in Section [5] using a sophisticated Markov Chain 
Monte Carlo technique. 



3 Gamma Rays from the Interaction of AGN Jets with the WIMP Halo 

In this section, we will compute the necessary components in order to predict the flux of gamma 
rays originating from the interaction of AGN jet particles with halo WIMPs. For the remainder 
of this paper, we will assume that the jet particles are exclusively electrons. The gamma ray flux 
originating from these types of interactions can be expressed as an integral over the energy of the 



impinging jet particles (Ee) as: 

dE, , (6) 



d<i>j f . (I d^^^^\ I 1 d^(Je+BH^-<+e+BH 

— OdM X 



di?^ 7 V'^yiGJV '^-^e / \ Mb dQdEj 

where the first term in the integral is related to the density of dark matter near the AGN core, 
the second factor accounts for the dynamics of the jet particles and the third factor is the cross 
section for the interaction of the jet particles with the halo WIMPs. Note that the first two factors 
(which depend both explicitly and implicitly on the characteristics of a particular AGN such as 
its distance from Earth, dAGN) determine the number of particles of each kind that can interact 
and the third factor (which depends on the scattering angle ^o between the jet axis and the line 
of sight) determines how frequently they interact at a given energy for the jet particle. In the 
subsequent sub-sections, we will consider each of the three above factors in detail. The discussion 
of the first two factors will follow closely the discussion in Rcfs. [2] and [3] and is only included here 
for completeness (and as a sanity check for the authors) . Those readers familiar with the results of 



Refs. [2] and [3] can safely skip to Section 3.3 



3.1 Dark Matter Density Profile 

The first factor in Eq. [6] is defined by the line of sight integral (as observed by the particles in the 
jet) of the dark matter density profile pDA-iir): 

5dm = / Pdm (r) dr , (7) 

where r„i„ is the minimum distance of interest (i.e., the "base" of the jet) and tq is the distance 
at which the jet becomes incohesive and irrelevant. Previous studies have found that, while results 
depend sensitively on the value of rmim the actual value of vq plays little or no role (since dark 
matter density profiles typically fall off steeply with increasing radius) . 

Although little is known about the dark matter density profiles of active galaxies, there are several 
models on the market. In this work, we focus on one particular model proposed by Gondolo and 
Silk [TW . In the Gondolo-Silk model (which assumes collisionless dark matter), the central black 
hole is assumed to grow adiabatically by the accretion of surrounding gas and stars. The result is a 
dark matter density profile with a very dense central spike. Assuming an initial dark matter density 
distribution with a power-law profile (p(r) ^ r"^), the Gondolo-Silk profile takes the form: 

The density in the core is dominated by Pcore which depends sensitively on the mass of dark matter 
and the age of the central black hole {tBu) as Pcore — Mum / ((o'w)o tsH)- The other component 
of the Gondolo-Silk profile is given by: 

p'W=Po(^)"(l-'^)7^)"'' (9) 



where 7sp is the slope of the density profile of the central spike, Rsp is the radius of the central 
spike and Rs is the Schwarzschild radius of the black hole. The spike slope and radius are given 
respectively by: 

7sp = -j^ ' ^-P = "7 ''o [^^) (10) 

and a-y oc 7''/^ [3j- Note that, because no stable orbit exists within the radius r < 'iRs, Pdm (r < 4 Rs) 
0. This condition sets a natural lower limit of integration Tmin = 4 Rs for Eq. [7] Finally, to fix the 
value of po ^ the normalization condition: 



/ ^Tir'^PDM {r) dr < AMbh , (H) 



is used. In other words, we require that the mass coming from dark matter be (at most) as large 
as the uncertainty in the mass of the central black hole (AMbh)- For the remainder of this paper, 
we will concentrate on the Centaurus A AGN. The parameter choices we use are [T7] : 



Mbh [ Mq ] Black Hole mass (5.5 ± 3.0) x 10^ 

Rs [pc] Schwarzschild radius 5 x 10^^ 

tBH [yr] Age of Black Hole 10* - 10^° 

aj 0.1 

To [kpc] upper limit of integration 15 

dAGN 3.7 Mpc 

6>o 68° 

The dark matter density profile for Centaurus A is shown on the left side of Figure [T] including 
four different scenarios depending on the value of the age of the black hole and the dark matter 
annihilation cross-section (ranging from 10~^° — 10~^^ cm^/s) for a WIMP with a mass of 100 
GeV. On the right side of the same figure, we also show the line of sight integral (Sdm) for the 
same four scenarios. As was first noted in [1], we see that a lower cross-section for pair annihilation 
allows a higher concentration of dark matter close to the spike. This enhances the probability for 
an interaction between a jet particle and a halo WIMP. We also see that younger black holes lead 
to a similar enhancement (due to there being less time for the black hole to deplete WIMPs living 
in the spike). Finally, we see that the eff'ective density that the jet probes varies very slowly with 
the lower limit of the integral (due to the spike in the density profile) , but decreases rather quickly 
outside of the spike region. 

3.2 The AGN Jet Factor 

Next, we discuss the second factor in Eq. ^ which accounts for the dynamics of the jet particles. 
As pointed out in Refs. [2] and [3], the exact geometry of the AGN jet is not very important for the 
purposes of our study. In fact, the most important factor affecting the behavior of the jet for the 
interaction studied in this work is the energy of the particles in the jet. Using the blob geometry 
as defined in Ref. |18] in which particles move isotropically and assuming that the gamma factor 
with respect to the black hole reference frame is not too large (F^ ^ 3) a broken power law energy 
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Figure 1: Dark matter density profile for Centaurus A for different values of Pcore where 7 = 1, 
Mbh = (5.5 ±3.0) X 10'^ Mq, i?s = 5 X 10*^ pc and tq = 15 x 10^ pc. Left dark matter density 
profile and right integral of line of sight of dark matter density profile. 



distribution is obtained from Fermi LAT observations [17j : 
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with si = 1.8, S2 = 3.5, 7hr = 4 X 10^, 7; 
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X 10^ , ^'max = 10 ' and 7' = E'/rUe. The constant 



fee must be found from the kinetic power of the jet (Le)'. 
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where the black hole reference frame is associated with unprimed quantities and the bulk reference 
frame with primed oneaM 

The final form of the second factor in Equation |6] is then: 
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where /io parameterizes the jet coUimation. In the following, we adopt a value of fio = 0.9 to ensure 
a highly-collimated jet. 

3.3 The Cross Section 

Finally, we consider the third factor in Eq. ^ (i.e., the differential cross section for the process 
Bh -|- e^ — > Bh + e^ + 1)- The Feynman diagrams which contribute to Bh + e^ — > Bh + e^ + j 
are depicted in Fig. I2I where the virtual particle is taken to be ii^t^'^) (the first KK excitation of 
the electron which we will simply denote as E). We write the total amplitude for this process as 



^ Lf, is taken as the Eddington limit for the black hole mass. 
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Figure 2: Feynman diagrams for the process Bu + e~ ^ Bh + e~ + 7 through £' exchange. 



■^ ~ X]i-^f^M(P5) where £^(^5) is the polarization vector for the photon. The amphtudes for the 
s-channel diagrams are given by: 
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where: 
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S12 - Me - ^^/sl2^£; , 



S34 - M% - i^/s34r_E , 



(17) 



(18) 
(19) 



and the kinematic invariants are defined by s.y = {pi +Pj)^ and tij — {pi ~pj)^. The energy- 
independent decay width F^; for E is computed to be: 



^E - 1^{Y^+Yr) j^ 



327r 



(20) 



The w-channel amphtudes are: 
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u{pi){YLPR + YRPL){h- h + Me){YlPl + YrPr){J,2~ h)rn(p2) , (21) 



u(p4)7''(i'4+ h){YLPR + YRPL){h- h + MeWlPl + YrPr)u{p2) , (22) 



M 



Uz 



^14^23 



1{va){YlPr + YrPl){^4- i>i + MEh^ih- h + Me){YlPl + YrPr)u{j>2) 



(23) 



The differential cross section for Bh{pi) + e{p2) -^ Bnips) + e{p4) + 1^{P5) is given in general 
by: 



da = {2ii)'^5'^(pi +P2-P3-Pi-P5) 



1 



(fp3 (Ppi 



d^P5 



4E1E2 {2Tr)^2E3 {2T:f2Ei (27r)32^5 ^^ 



\M[^ 



(24) 



Using the 5'^{pi-\-p2—p?,—PA~Ph) to integrate over d^p^ and the remaining 5{Ei+ E2 — E^— Ei^ — E^,) 
to integrate over dE^, we arrive at a differential cross section of the form: 
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(25) 



We choose to work in the initial WIMP rest frame and align our coordinate system such that the 
initial electron four-momentum lies along the -|-z-axis and the azimuthal angle for the photon is 
zero (05 = 0). In this frame, the relevant momenta are given by: 



pi = {Mb,0) 

P2 = ^2(1, 0,0,1) 

P4 = -£4(1, sin 04 cos (/)4, sin 6*4 sin 04, COS 514) 

P5 = ^5(l,sin6'5,0,cos6i5) 



(26) 
(27) 
(28) 
(29) 



and we use conservation of four-momentum to replace p^ (i.e., p^ — pi + P2 — Pi ^ Pb)- Note that, 
for Cen A, 9^ = 9q = 68°. The energies for the final state WIMP {E^) and electron (£^4) in this 
frame are given respectively by: 
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Mb -I- i?2 (1 — cos ^4) — £'5 (1 — sin 9^ sin 94 cos 04 — cos 9^ cos ^4) 
while the Jacobian J takes the form: 

J — —- [E4 — E2 cos 04 + £5 (sin 05 sin 04 cos ^4 -f cos ^5 cos 6*4)] . 



(30) 
(31) 

(32) 



The differential cross section as a function of the incoming electron's energy [E^ — E2) \s shown 
in Fig. [3] for one particular photon energy {E~f = E^ — IQ GeV). For this point in phase space, 
we obtain two resonances: the first (located near Me — Mb — 15 GeV) is due to pieces of the 
amplitude-squared which go like ~ I/IE12P, while the second resonance originates from pieces of 



Differential Cross-Section 




£,= 10GeV 
Mb=200 GeV 
Mi=215 GeV 



50 70 

£, (GeV) 



Figure 3: The cross section for Bh + e~ — >■ Bh + e~ + 7 for photon energy E^ = 10 GeV as a 
function of the incoming electron energy. The resonance located near 15 GeV is due to pieces of 
the amplitude-squared which get contributions from diagrams ^2 and ^3 (which depend on the 
propagator term E12), while the second resonance located near 27 GeV originates from diagrams 
5*1 and S^ (which depend on the propagator term E34). 



the amplitude-squared which go like 
over the second. 



1/1^34^. We note that the first resonance always dominates 



We are now in a position to compute the gamma ray flux from interactions of AGN jet particles 
with halo WIMPs. In the following, we will actually present our results in terms of the spectral 
energy distribution, vS^, which is simply related to the flux as: 



vS^ 



p2 X rf*7 



(33) 



Using Eqs. Q, ( [l4| ) and (25) in Eq. ^ and integrating over the initial electon's energy, we obtain 
spectra such as that shown in Fig. |4] For comparison, we separately calculated the results for the 
5-d UED case which were previously done in Ref. [2]. We note that, in both cases, the spectrum 
rises sharply as a function of the photon energy until a sharp cutoff at the mass difference between 
the WIMP and the up-scattered particle (i.e., the E particle in the 6-d case). The sharpness of 
this cutoff will be softened slightly once detector resolution effects are taken into account (see 
Section [5]). 

Finally, we should note that we have performed this calculation in two ways. In addition to 
calculating the photon flux using the exact kinematics as outlined above, we also computed the flux 
using the collinear approximations proposed in Ref. [2] . We found excellent agreement between the 
two calculations in the regions of the resonances and only small differences outside these regions. 
Since the resonances give the largest contributions to the flux, we chose to use the more compact 
collinear- approximated calculation for our scans described in Section [5] since these resulted in much 
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Figure 4: The spectral energy distribution vS^, (= E^ x d'^^/dE^) versus the photon energy from the 
interaction of WIMPs with the AGN jets of Centaurus A. Note that 5D case refers to a comparison 
to the model studied in [2] and 6D case corresponds to the model studied in this work. 

faster computations. The details of the cross section calculation in the coUinear approximation are 
summarized in Appendix [A] 



4 Gamma Rays from Dark Matter Annihilation 

In this section, we compute the flux of gamma rays coming from WIMP annihilation. In contrast to 
the flux from jet-halo WIMP interactions, the flux from WIMP annihilations depends quadratically 
on the density profile (because two WIMP particles are involved in the initial state) and is given 
by: 



d$ 



dE^ 
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dE^ SnM^d^Qj^ j^^ 



where {av)tot is the total annihilation cross section and dN^/dEj is: 
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(34) 



(35) 



We use the index / to denote the annihilation channels with one or more photons in the final 
state, {<Jv) f is the corresponding cross section and dN^ /dE^ is the normalized photon spectrum 
per annihilation. 

Below, we briefly outline the calculation of the two contributions to the WIMP annihilation gamma 
ray spectrum: (z) the continuum from annihilations into SM particles which then subsequently 
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Figure 5: Gamma ray spectra from continuum WIMP annihilations. Left: The differential photon 
spectrum per annihilation event as a function of the fractional photon energy x = E-^/Mb for 
Mb = 200 GeV. Right: The spectral energy distribution vSi, (= E'i^ x d^j/dEj) versus the photon 
energy for Centaurus A for two choices of the dark matter density profile. 

radiate or hadronize/decay into photons and (ii) line emissions from loop- induced decays directly 
into photons. For more details, the interested reader should refer to Ref 6^ where these calculations 
were first performed. 

4.1 Continuum 

As discussed earlier, pairs of LKPs annihilate predominantly into pairs of electroweak bosons WW 
and ZZ, SM Higgs bosons HH and (if heavy enough) top quark pairs tt. We compute the continuum 
gamma-ray spectrum from these annihilations using the micrOMEGAs code [15]. For Mh = 
125 GeV and Mb = 200 GeV, the annihilation fractions are roughly 48% BhBh -)■ WW, 28% 
BhBh -^ HH, 22% BhBh -^ ZZ and 2% BhBh -^ ti and the total cross section is {av) = 
2.40 X 10-26 cmVs. 

In Fig. [5] (left), we plot the differential photon spectrum dN/dx where x = E^/Mb- In Ref. |5], it 
was pointed out that, in contrast to other models, the continuum spectrum from the chiral square 
model sharply decreases well before the value of the WIMP mass Mb ■ The reason is that the main 
annihilation modes are dominated by "photon unfriendly" modes, consisting of massive (and often 
neutral) particles which are unlikely to radiate high-energy photons. The result is that the bulk 
of photons come from radiation (or after hadronization, decays of vr^'s) from the even softer decay 
products of the particles produced in the primary annihilation. 

Fig. [5] (right) depicts the spectral energy distribution for continuum gamma rays from WIMP anni- 
hilations for Centaurus A using two different dark matter density profiles (the Gondolo-Silk profile 
discussed above and the Gnedin-Primack profile of Ref. [10] )■ Comparing to the spectrum from 
jet-halo interactions, we see that below ~ 15 GeV the jet-halo gamma rays dominate by an order 
of magnitude, but above this value of E^ the continuum from WIMP annihilations is of the same 
order as the jet-halo contribution. Fig. [6] shows the sum of the jet-halo and continuum annihilation 
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Figure 6: Total spectral energy distribution for the process e + B with one or two final photons 
(leading order and one loop) including continuum emission from dark matter annihilation. 
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Figure 7: Example Feynman diagrams for the process Bh + Bh -^ ^ -^ V . Q 1^^ are the gauge 
eigenstates for the (1,0) excitations of the SM leptons {£). 



gamma ray spectra for a particular choice of Mb and the mass splitting {Me — Mb)- 



4.2 Line Emission 

Let us now consider the direct (loop-level) annihilations into gamma rays, 7 + X. In the chiral 
square model, the WIMP is a scalar particle such that the only possibility is for X to be a vector 
particle. This allows for the possibility of three lines originating from the 77, 'yZ and 7!?'^^'^) final 
states (where B^^'^^ is the (1,1) KK excitation of the hypercharge gauge boson with mass roughly 
~ \/2Mb)- In fact, the chiral square model has been shown to produce a gamma ray spectrum 
from our galactic center that possesses the rare characteristic of multiple and well-separated line 
features (even after the inclusion of detector resolution effects) [6] . 

Sample Feynman diagrams which contribute to Bb(pi)Bh{p2) -^ l{PA)y{Pv) are depicted in Fig. 



13 





10 " 


.,...,..., 











<"l, ^> 


tn 


TQ-28 






E 


10-^^ 


^\r-~.l_^~~~~~~' 


._^ _ 1 


A 
> 


TO-30 


- '■■■'■...^ ^ ~--- 


■— _j 


b 
V 


10-^' 




" ""■■■" 




10-^^ 




1 



200 



Figure 8: The cross sections for BhBh 
This plot was taken from Ref. [6]. 
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7Z and 7i3(^'^) as a function of the WIMP mass. 



[7I The amphtude for this process can be written as: 

M = e^^*ipA)e'gipB)M^"'ipi,P2,PA,PB) 



(36) 



where e^ and e^ are the polarization tensors of the photon and V gauge boson, respectively. In 
general, the tensor A^^"^ can be expanded in terms of the metric tensor and external momenta 



M^ 



Ag^'' + B.p^^p'i + B^pt^p^, + Bsp'ip^, + B^p^^pl + B^p^p^^^ 
+Beptp\ + Syp^p^ + Bsp^^p\ + B^p^i^pl . 



(37) 



It was shown in Ref. [5], that using the non-relativistic nature of the WIMPs (i.e., pi ~ p2 = 
{Mb,0)) along with conservation of momentum, the only surviving term is g'^'^ and, thus, the 
computation of the cross section is very much simplified. Fig. [8] shows the cross sections for 77, 7Z 
and 75^^'^-' as a function of Mb for the KK masses AIe = 1.2AIb and A/^d.i) = 1.6Mb- 

The enhancement of the 75'"'^'"'^^ cross section with respect to the other two may seem surprising 
due to the fact that the B^^'^^ is so much more massive than the photon or Z. However, the effect is 
understood as follows. In the case where V is either a photon or a Z boson, the couplings between 
Vii and V^^ (where ^ are the gauge eigenstates of the (1,0) excitations of the SM leptons) are nearly 
the same strength and the result is a significant cancellation between the various diagrams. However, 
in the case where V is identified with the B^^'^^^ KK gauge boson, the Vil couplings are loop- 
suppressed, while the V£,^ couplings are relatively large. This results in less cancellation between 
the various diagrams and an enhanced cross section compared to the other two processes. 

The spectra for the lines will depend on the intrinsic width of the particle produced in association 
with the photon. In other words, while the 77 line will simply be a delta function at the WIMP 
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mass Mb, lines arising from the process BhBh — > "fV will possess an intrinsic width which will 
depend on the mass of the boson in the final state My : 



where Ty is the width of V and: 

A ^ 

/2 = 



dN^ _ AMbMvVv 



tan 



Mv 
Mf_ 



1 /4:Ml~M^ 

tan^i ^ ^ 

V MvTv 



{AMI ' ^MbE^ - Mlf + r^M^ 



(38) 

(39) 
(40) 



In practice, the detector resolution tends to dominate the line shape for the cases we consider. 



5 Results 

In this section, we combine the fluxes from jet-halo WIMP interactions with those of WIMP anni- 
hilations (both continuum and line emissions) and compare with the data collected by the Fermi 
gamma ray telescope [17'. In addition to these fluxes, we add a background contribution assumed 
to originate from astrophysical processes that behaves as a power law: 



d$ 
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'^{ih) 



(41) 



We assume the 7-ray spectra is unaltered in route to measurement and we smear the energy ac- 
cording to a Gaussian kernel: 



d$t 



dE^ 



G{){E E^) 



fd<^ 



bkg 



d$. 



V dE' 



dE' 






AGN 



dE' 



dE' 



where the kernel is: 
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2Tra, 



(42) 



(43) 



exp 



Gexp is the experimental resolution, which we adopt as a exp — 0.1-E' for Fermi LAT. The scale 
factor Aagn parameterizes the uncertainty in the AGN contribution, which we take to be 0.5 — 1. 
To arrive at the x^ contribution from the AGN data, we average the event rate over the energy 
bins. 

To fit the model, we perform a Bayesian fit to the data, the details of which are provided in Ap- 
pendix Ib] The fit includes the Centaurus-A AGN data PTj. We fit the dark matter relic abundance 
in one of two ways: 

^B«/i^ = 0.11 ±0.01, (44) 

which we denote as the "saturated" case, where the relic abundance is explained entirely by the dark 
matter candidate. Here, we take a 10% uncertainty that is prescribed to be the theory uncertainty 
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on top of the excellent experimental measurements of Planck ,21j. In the other case, we consider 
only the measurement as an upper bound of the dark matter relic abundance which we refer to as 
the "unsaturated" case. We will see that these two circumstances lead to different characteristic 
spectra, and therefore a different overall fit to the AGN data. 

To further constrain the model, we will also include in our fit the most current limits on WIMP 
parameter space from direct detection experiments (in particular, those from XenonlOO |22j). Since 
we allow for a reduced relic abundance, we scale the scattering rate exclusion and Bh annihilation 
observables with the relic abundance: 



_incas. 



nR„h^ 



fsi 



fi. 



h-^' 



^x 



av^x 



nR„h^ 



fin 



h^ 



(45) 



where the measured rate must be compatible with present experimental bounds from Fermi LAT ^ 
and XenonlOO |22l. 
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Figure 9: Photon spectra from the best fit point for the saturated (left) and unsaturated (right) 
rehc abundance scenarios. The unsaturated rehc abundance case yields an energy cutofi" that fits 
the Centaurus-A data slightly better than the saturated case. While the continuum and lines are 
generally suppressed below the astrophysical background, the B^^'^^^ line can appear well above 
the continuum, off"ering possible detection in our own galaxy. 

In Fig. [9) we show the 7-ray spectra from the best-fit of the Centaurus-A data for the two relic 
abundance scenarios we consider. For the astrophysical background, we find the power law parame- 
ters for the saturated case to be Ab — (6.1^Qg) x lO"-"^^ and 5h = — 3.06=b0.14, with the unsaturated 
scenario having a very similar fit. We find that the saturated relic abundance case yields a cutoff 
that is of higher E~^ than in the unsaturated case. We observe that since this cutoff is related to 
the mass difference 

£;-t=M^-AfB, (46) 
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the model is forced into a particular and well defined relic density scenario which we will discuss 
below. Generally, the continuum from dark matter annihilation and the 7-ray lines are suppressed 
well below the astrophysical background. However, the B'^^'^'>^ line can appear well above the 
continuum, allowing detection within our own galaxy. 
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Figure 10: Dark matter relic abundance and mass posterior distributions while requiring the DM 
candidate is not above the measured abundance (top left panel), and fully saturates the measured 
abundance (top right panel). The bottom panels show the relic abundance and mass splitting 
illustrating the role of Bh — E co-annihilation. The dashed (dotted) contour contains the la (95%) 
C.L. regions, respectively. 



In the top panels of Fig. 10 we show 2-d posterior distribution of the relic density and dark matter 
mass under the two above scenarios. In the unsaturated case, there is a preferred region at nearly 
half that of the measured relic abundance with mass in the 180 — 260 GeV range. However, in the 
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other case, the rehc density is indeed saturated, but with a lower dark matter mass of 160 — 200 
GeV. The source of the lower relic abundance in the first case is due to coannihilations of Bh with 
the compressed spectrum of the first excitations. This is highlighted in bottom panels of Fig. [Toj 
which shows the relic abundance and mass splitting of the excited electron. As previously seen, 
this is important as the "edge" in the Fermi data point to a mass splitting on the order of 15 
GeV. A spectrum with a mass splitting this small has a significant coannihilation rate, thereby 
lowering the relic abundance. Alternatively, a non-standard cosmological scenario such as non- 
thermal mechanisms that enhance the relic abundance must be realized. 
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Figure 11: Scattering cross section for dark matter mass in the unsaturated (left) and saturated 
(right) relic abundance cases. The dashed (dotted) contour contains the Itr (95%) C.L. regions, 
respectively. Either scenario is easily probed by the LUX and Xenon IT experiments. 



Finally, in Figs. 11 and 12 we present associated dark matter observables. In Fig. [TTj we see 
the spin-indpendent scattering cross section for the two relic abundance cases. In the case of a 
saturated relic abundance, the expected cross sections are quite close to the XenonlOO limit [22]. 
This indicates some amount of tension within the model. However, the unsaturated case is naturally 
within a factor of 3-4 of the current bound. Either case is within reach of future direct detection 
experiments such as LUX [M" and Xenon-lTon ^^5] with a cross section no lower than 0(5 x 
IQ-io pb). 



In Fig. 12 we show the location and predicted cross section into the 7-ray lines 77, Z7 and B^^'^'j 
and the present 4 year Fermi line search limit |53]. The B'^^-'^'>^ is sizable and is well positioned for 
Fermi line searches within our own galaxy in the near future. We find the B^^'^'j should have a 
line location of 50 — 100 GeV and cross section at the 0(10^^* cm^/s) level. The 77 and Z7 lines 
are of similar strength. Furthermore, given the energy resolution of 7-ray experiments, they will 
likely be smeared together as one line at this energy. Therefore, they are combined in the contour 
plot as one line. 
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Figure 12: Cross section distribution of 7-line from BuBu -^ X + ^ with and without the relic 
density saturated. The dashed (dotted) contour contains the la (95%) C.L. regions, respectively. 
The line cross section limit from Fermi after 4 years is given by the dashed line above the contours. 

6 Discussion 

Given the results of the fit to the Fermi data for Centaurus-A, we now discuss further implications 
consistent with these data. Since the best fit points in the model correspond to a mass splitting of 
0(15 GeV), we anticipate that the relic density of Bh should not necessarily be the value measured 
by WMAP and Planck due to coannihilation with the up-scattered state, E. Therefore, a second 
dark matter candidate must be needed to explain the remaining dark matter of the universe; a simple 
explanation is the axion arising from the Peccei-Quinn solution of the strong CP problem. 

With a small mass splitting between the dark matter candidate and the compressed spectrum 
associated with the first level excitation, we expect searches at the LHC to have difficulty discovering 
this scenario. This is due to the soft decay products that fall below the cuts required in such analyses. 
Indeed, the null searches for new heavy states in models such as supersymmetry have not precluded 
this scenario. However, a mass scale of 0(200 GeV) is within reach of a future ILC at 250 or 500 
GeV. There, the cleaner background and more precise control over initial state energy may allow a 
discovery. 

It is worth noting that in the presence of brane kinetic terms, the masses in extra dimension models 
are allowed to vary more freely. In this case, while the strength of coannihilation between B^j 
and the compressed spectrum will be weakened, we can say that coannihilation between Bh and 
E will take place to some level due to the AGN jet data requiring a 10 — 20 GeV mass splitting. 
However, the searches at the LHC may uncover the decay of, say, Q^^^^^ -^ Bnq with a sufficient 
mass splitting. 

Future data from many experiments should further probe this scenario. The Fermi LAT should 
be able to shortly discover the B^^'^'^ line. Moreover, future direct detection experiments such 
as LUX [21] and Xenon- ITon [55] are well positioned to probe the entire scattering cross section, 
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assuming the Bh relic abundance is not suppressed too much. The Fermi LAT and HESS mea- 
surements on Centaurus-A should also provide additional statistical power to confirm or refute this 
scenario directly. 



7 Conclusion 

With the range of viable dark matter parameter space continually shrinking and the validity of 
certain "anomalies" still under debate, it is crucial that the search for dark matter be extended 
to ever more exotic scenarios and/or signal sources. One such possible source may be the core of 
active galactic nuclei where the density of dark matter is expected to be extremely high. Recently, 
several groups have explored the possibility of observing signals of dark matter from its interactions 
with the high-energy jets emanating from these galaxies with favorable results. 

In this work, we built upon these analyses by including the other components of the WIMP-induced 
gamma ray spectrum of active galactic nuclei; namely, (i) the continuum from WIMP annihilation 
into light standard model states which subsequently radiate and/or decay into photons and [ii) the 
direct (loop-induced) decay into photons. In particular, we considered the gamma ray flux from 
the so-called "chiral square" model which is a model of two universal extra dimensions. This model 
was chosen because of the interesting annihilation features it possesses (i.e., a "forest" of gamma 
ray lines). 

We found that, like other models considered previously, the WIMP-induced gamma ray flux of 
AGN from the chiral square model exhibits several interesting features which may be observable 
with the Fermi gamma ray telescope. In addtion to this, we have shown that, in conjunction with 
other measurements from direct and indirect detection experiments, the gamma ray spectrum from 
AGN can be an invaluable tool for restricting WIMP parameter space. 
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A Cross Section in the Collinear Approximation 

As is evident from Fig. [3] the dominant contributions to the cross section and, thus, the gamma ray 
flux come from the s-channel Feynman diagrams due to resonances when the intermediate E goes 
on-shell. Thus, it appears a safe assumption would be to neglect the u-channel diagrams shown in 
Fig. [2] since there is no chance for the intermediate E in these diagrams to go on-shell. 

Besides the resonant-enhancement from the E going on-shell, there is another source of enhancement 
in the process considered here: a logarithmic enhancement when the photon is collinear with the 
final-state electron (note that, since the angle of observation does not line up with the jet axis, there 
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is no possibility of the photon being colhnear with the initial-state electron). Below, we outline 
the calculation assuming that the resonant and coUinear regime makes up the most important 
contribution to the cross section. This calculation and this appendix follow closely Appendix B 
of Ref. P] (by design) and provides an important check on our calculation using the full set of 
Feynman diagrams and the exact kinematics of the process. 
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We group the pieces of the amplitude-squared that contain the collinear logarithm (|A4|fo„) and 
those that do not i\M\'^gi ) as: 
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B Bayesian fit using Markov Chain Monte Carlo 

In this work we utilize a Markov Chain Monte Carlo (MCMC) to fit the available data. The MCMC 
approach is based on Bayesian methods to scan over specified input parameters given constraints 
on an output set. In Bayes' rule, the posterior probability of the model parameters, 0, given the 
data, d, and model, Af , is given by 

where 7r(0|M) is known as the prior on the model parameters which contains information on the 
parameters before unveiling the data. The p{d\9,M) term is the likelihood and is given below in 
Eq.l57J The p{d\M) term is called the evidence, but is often ignored as the probabilities are properly 
normalized to sum to unity. A common algorithm for MCMCs is that of Metropolis and Hastings 
(MH), in which a random point, 9i, is chosen in a model's parameter space and has an associated 
likelihood, £,;, based on the applied constraints. A collection of these points, {9i}, constructs the 
chain. The probability of choosing another point that is different than the current one is given 
by the ratio of their respective likelihoods: min(-^,l). Therefore, the next proposed point is 
chosen if the likelihood of the next point is higher than the current. Otherwise, the current point 
is repeated in the chain. We adopt the likelihood 

Q ^ e-^.x]!i ^ g-s,te,-d,)V2^|^ (57) 

where yij are the observables calculated from the input parameters of the i*'' chain, dj are the values 
of the experimental and theoretical constraints and Oj are the associated uncertainties. 

The advantage of a MCMC approach is that in the limit of large chain length the distribution of 
points, Qi^ approaches the posterior distribution of the modeling parameters given the constraining 
data. In addition, the set formed by a function of the points in the chain, f{9i), also follows the 
posterior distribution of that function of the parameters given the data. 

Instead of the MH algorithm, we adopt the Goodman- Weare (GW) algorithm [35], which has been 
shown to be more efficient at choosing trial points and therefore increases speed. A collection of 
initially random points referred to as "walkers" at time t, {9i(t)}, are placed in parameter space. 
A given walker's next trial position, 9i{ty, is based chosen along a line connecting it to another 
randomly chosen walker, 9j, with the trial 

9,{t) ^ 9'^ = 9, + Z i9,{t) ^ 9,) , (58) 

where Z is a random variable with the distribution 

giz) (X ^ (59) 



in the interval [l/a,a] and vanishes elsewhere. The trial point is accepted if 

q = z —- — > r, (60) 



where r is a uniform random variate. This selection ensures detailed balance within the algorithm, 
which is necessary for the algorithm to be statistically relevant. The parameter a sets the scale 
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of jumps. The value a = 2 is often used as it yields a very efficient exploration of the parameter 
space. 

We allow a "burn-in" period to prevent correlation of the walker's initial random state to the set of 
walkers used for the posterior distribution. We guarantee this by constructing the autocorrelation 
time 

am = — - Yl (^'(*) - ^'*)(^^(^ + ^) - ^'^) (61) 

for each input parameter and verify that the burn-in length is at least 5 times the auto-correlation 
scale. 
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